Entanglement and area law with a fractal boundary in a topologically ordered phase 
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Quantum systems with short range interactions are known to respect an area law for the entanglement entropy: 
the von Neumann entropy S associated to a bipartition scales with the boundary p between the two parts. Here 
we study the case in which the boundary is a fractal. We consider the topologically ordered phase of the toric 
code with a magnetic field. When the field vanishes it is possible to analytically compute the entanglement 
entropy for both regular and fractal bipartitions {A, B) of the system, and this yields an upper bound for the 
entire topological phase. When the A-B boundary is regular we have S/p — 1 for large p. When the boundary 
is a fractal of Hausdorff dimension D, we show that the entanglement between the two parts scales as S/p = 
7 < 1 /D, and 7 depends on the fractal considered. 

PACS numbers: 03.65.Ud, 03.67.Mn, 05.50.+q 



Introduction. — Entanglement is certainly one of the most 
striking aspects of quantum theory. Not only is it the key in- 
gredient for protocols ranging from quantum teleportation to 
cryptography, but it also has an important role in the study 
of condensed matter and many body systems Quantum 
phase transitions can be understood in terms of entanglement 
P], and new exotic states of matter that defy a description in 
terms of local order parameters show a signature of topolog- 
ical order in the global pattern of their entanglement 
Moreover, the analysis of the scaling of entanglement in the 
ground state of condensed matter systems has shed new light 
on the question of their simulability [5]. 

Especially for the last reason, one is interested in knowing 
how entanglement scales with the size of the system. If there 
is a gap, all correlations decay exponentially with the distance 
in units of the length scale [j6fl. In this case, one also expects 
the entanglement to be short ranged, so that only the degrees 
of freedom of the boundary of the system contribute to the 
total entanglement. This is the so called area law for the en- 
tanglement (see Ref. jltl for a comprehensive review). 

In this work, we study the case of a topologically ordered 
state, the ground state of the toric code ||9|]. For this state - 
and a class of topologically ordered states - the entanglement 
can be computed exactly 1^. For a bipartition with a regular 
boundary p, the entanglement measured by the von Neumann 
entropy S is exactly S — p — 1, where the coiTection —1 is 
due to a topological contribution to the entanglement 
Obviously, 7 := S/p is 1 in the limit of large p. If we add 
perturbations to the model, topological order is not destroyed 
until a quantum phase transition happens. Throughout the en- 
tire topological phase the entanglement is upper-bounded by 

t ' 

its value in the unperturbed model [10]. 

Here we study the case in which the boundary of the sys- 
tem is a fractal curve of Hausdorff dimension D. This situa- 
tion arises under a large variety of experimental conditions in 
two-dimensional systems flill . The scaling of entanglement 



for self similar systems is important also in view of devis- 
ing efficient algorithms which use the renormalization group 
for computing ground states of quantum systems in two di- 
mensions [5]. One could expect that as the boundary of the 
system becomes less regular, the entanglement increases with 
the length p of the boundary, as in the case of fermions lfl9ll . 
In contrast to the fermionic case, we find that for topologi- 
cally ordered spin systems the entanglement decreases with p. 
The length of a fractal curve - and consequently the entan- 
glement - diverges in the limit of exact fractality [13]. How- 
ever, for every step n of the iteration of the fractal, the length 
of the curve is a finite number which increases with n. 
In contrast to regular boundaries, for fractal boundaries 7 is 
a fractional number: we can speak of fractal entanglement. 
Moreover, we shall see that 7 < D^^. 

Entanglement and topological order. — Consider a unitary 
representation of a group G acting on spin- 1/2 degrees of 
freedom with Hilbert space Ti.. Since we wish to compute the 
entanglement entropy associated to a bipartition of the system, 
we are interested in the properties of the group when we split 
the Hilbert space as = T-La ® T-Lb- We assume that there 
exists a product state |0) = \Qa) ® |0s) G Tf.. We can now 
define the (normalized) G-state as |^g) '■— J2geG '^(s)^!^)- 
If all the coefficients are equal, we call the state a G-uniform 
state: \G) := \G\-^/^ J^g^G s\^)' where |G| is the order of 
G. Note that \G) is stabilized by the group G. Let us now 
define the two subgroups of G that act trivially on the subsys- 
tems A, B respectively: Ga '■= {g & G \ 5 = ® Is} and 
similarly for Gb- By defining the quotient group Gab '■= 
G/(Ga X Gb), we can write G as the union over all ele- 
ments of GAB^■ G = \J[h]eGAB^^9A ^ gB)h \ gA ® T^b e 
Ga, ^a ® 9b £ Gb}- The state can thus be written as 
*g) = \G\^^^^ J2gA^gB&GAxGB a{gA®gB-,h)hA®hB® 

heGAB 

{9A ® 9b)\0)- If the coefficients a in the expression for 
\'^g) satisfy the separability condition a{gA ® gB,h) = 
a{gA ® gBhg) = aA{gA)aB{gB)l3{hg) for every g e G, 
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then it is possible to prove fl^ that the von Neumann en- 
tropy of the G-state corresponding to the bipartition {A, B) is: 
S{\^g)) - -Y^MeG^u \NANBm?^^g2\NANBm\\ 
where Nj, := Eg^eG^ Mgx)?. for X = A,B. By con- 
vexity of 5 we have5(|*G)) < S{\G)) ^ log^ \Gab\- 

This formaHsm is remarkably well suited to describing 
topologically ordered states. In many quantum spin systems, 
topological order arises from a mechanism of closed string 
condensation and the group G is the group of closed strings 
on a lattice llljll . An important example of topologically or- 
dered system is given by Kitaev's toric code, which provides 
a model for which at zero temperature topological memory 
and topological quantum computation are robust against arbi- 
trary local perturbations [9]. The model is defined on a square 
lattice with spin- 1/2 degrees of freedom on the edges and pe- 
riodic boundary conditions. To every plaquette p we associate 
the operator product of cr^ on all the spins that comprise the 
boundary of p, i.e., Xp — Wjfzp^^^j- To every vertex s we 
associate the product of cr^ on all the spins connected to s: 
Zg = rijes "^i ■ operators Xp generate a group G of 
closed string-nets. The Hamiltonian of the toric code in an 
external magnetic field is: 



Hu 



= -Y,Xp-\Y,Zs + {i-\) 



(1) 



where we have introduced a control parameter A. A second 
order quantum phase transition at Ac ^ 0.7 separates a spin- 
polarized phase (0 < A < Ac) from a topologically ordered 
phase (Ac < A < 1) |1^ 15]. The ground state of iJtoric 
is a G-state throughout the entire topological phase. It is G- 
uniform at the toric code point A = 1, and becomes less uni- 
form as A decreases to Ac. 

We now wish to argue that the separability condition for 
a{g) is satisfied throughout the entire topological phase, and 
hence by convexity 5a < S{\G)) — logj \Gab\ for Ac < 
A < 1, with the bound saturated at the toric code point. At 
A = the ground state is the uniform superposition of closed 
strings. The A term in Eq. ([TJ is a tension for the strings. 
As we increase A, larger strings become less favored in the 
ground state. Everywhere in the topological phase, that is, 
for sufficiently small A, the ground state is still the superposi- 
tion (with positive coefficients lUB]) of closed strings g G G. 
The expectation value (g) of any closed string g G G of 
length I (a Wilson loop) can be written as {g) = Gj^e'^^'^'''^^ 
where Gt is a constant that does not depend on g (due to 
translational invariance). Similarly, in the polarized phase we 
have (g) = Gpe^^°-''^\ where a is the area enclosed by the 
string [11]. Now, we know that {g) — \a{g)\'^ at any point 
in the topological phase, since the ground state is a G-state 
and does not contain any open strings. Since the length I 
for a given string g = gA ® gshg can be decomposed as a 
sum of the corresponding substrings, / = Ia + Ib + Iab, 



we have a{g) 



Gte-'(9)/2 = Gte 



Ia/2 -Ib/2 -Iab/2 ^ 



aA{gA)o:B{gB)P{hg), i.e., we have separability. 

Fractal boundary. — Henceforth we consider the toric code 
point A = 1, where S = log2 IGabI. We define biparti- 



tions by drawing strings along the edges of the lattice. One 
can prove |3] that log2 |Gab| is the number of independent 
plaquette operators Ap acting on both subsystems A and B, 
which in turn is the number of squares that have at least 
one side adjacent to the boundary p of the region A, see 
Fig. [T] How do we measure pi We shall show that the 
support of the mixed part of the reduced density matrix is 
given exclusively by the spins on the boundary. This mixed 
part is the only part contributing to the entanglement between 
the A and B partitions. Therefore we define the length p 
as the number of boundary spins. Indeed, letting Qx = 
\Gx\^^^^ 'Yg^^Gx 9^' "^^^^ X — A,B, the ground state can 
be written as |G) = \GAB\~^^^J2heGAB '^^'^bQaQbIO). 
It follows from the definition of Gab that we can pick Ha 
up to local transformations of the loops inside A and B. 
Specifically, we can pick Ha as acting only on the spins on 
the boundary. Since Qa,Qb are local operators, the re- 
duced density matrix of the ^-subsystem is equivalent to 
one separable as TrB[|G)(G|] — \ip){tp\ Pa, where 
is a pure state describing A's bulk, while the mixed part is 
PA = \Gab\^^ Y^heGAB hA\0){0\hA, where/iA acts exclu- 
sively on the spins along the boundary of A fllTll . Thus S/p is 
the average entanglement per spin in the support of pA ■ 

We now consider the case of a bipartition defined by a 
closed fractal curve. Since the model studied here is defined 
on a square lattice, we consider bounded regions of 7? de- 
pending on a parameter n, denoted by An- Here n repre- 
sents the number of steps in the iteration generating the fractal 
curve. The perimeter of An is denoted by p [An)- The num- 
ber of squares of size one adjacent to the boundary of An is the 
entanglement S{An) associated to the bipartition {An^Bn)- 
We are interested in the large n limit of the ratio between en- 
tanglement and perimeter: 7 {A) := lim„_j.oo S{An)/p{An)- 
One might expect the scaling law 5 = p — 1 to be independent 
of the geometric properties of the bipartition, but this is not the 
case. From Fig. [T] we see that when the boundary of A has 
some inward angles, or wells, or other "kinks", the number of 
squares adjacent to it is less than the length of the boundary 
around it. For instance, an inward angle, a well, and a hole all 
have just one adjacent square of side 1 but they have lengths 
2. 3, 4 in the lattice spacing unit, respectively. We call a and h 
the number of inward angles and holes, respectively. It is not 
hard to show that |jl8|] 



S ^ p — a — 3h. 



(2) 



We wish to study how these numbers scale for a fractal expan- 
sion, and find the corresponding scaling of the entanglement. 

In the following, we shall compute 7 for several fractal 
curves. The results are summarized in Table|l] The main result 
is that, depending on the fractal region, 7 can be a fractional 
number The Hausdorff dimension D of the fractal does not 
uniquely determine the value of 7, but (in all the examples 
considered) we have the bound 7 < D^^. 

Examples. — The Sierpinski carpet on 1? , denoted by iS„, 
is a bounded region of 1? defined iteratively in the following 
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FIG. 1 : The drawings show different bipartitions of the system. The 
subsystem A consists of all the spins marked by the black squares. 
The entanglement is given by the number of plaquette operators act- 
ing on both subsystems, marked by red dots. For a regular figure 
(left), this number coincides with the perimeter p, which is the num- 
ber of spins along the boundary (in yellow). Every time there is an 
inward angle, there is one such operator for three units of length. 
The well (middle) contains two inward angles. A hole (right) of size 
1 accounts for 4 units of length and contains only one star operator. 




FIG. 2: (Color online) Top, left to right: Sierpinski carpet <S3, Greek 
cross C/3, Minkowski sausage I3, T-square £4,- Bottom, left to right: 
Moore polygons M3, Vicsek fractal V3, half perimeter of the Koch 
polygon AC5, 4 X 4 chessboard C4. 



way: f/j iSi is a 3 x 3 square without the central 1x1 square. 
The Sierpinski carpet Si has a single square hole, (ii) Sn+i is 
a bounded region inscribed on a 3" x 3" square on Z^. This 
is obtained by placing 8 copies of Sn on all quadrants of the 
square, but the central one (see Fig. |2]i. Given the recursive 
structure of iS„, direct calculations show that a (iS„) — ]^8" — 
|. The number of equal holes of side 3* is 8"^^^', so h{n) = 
8"^^. Observe that the external perimeter of iS„ is 4 x 3". 
Then the perimeter is p(5„) = 4(3"+3"-i+X;r=o^(3' x 
8""!"*)) = 4 (4 x 3" + 8") /5. With this information, from 
Eq. ^ we obtain 7 (5„) = 99/224. 

The Greek cross on 1?, denoted by Qn, is a polygon in 1? 
defined by a closed path of length p{Qn) 8n + 8, including 
the point (0, n) and the step {(0, n), (1, n)}. The path max- 
imizes the number of inward angles over all the closed paths 
of the same length including the point (0,n). Fig. Ogives 
the first few instances. It is then evident that a(Qn) = An. 
For this polygon, h{n) = and thus from Eq. (|2|i we have 
S{n) ~ p{n) — a{n). Therefore, 7 {On) = 1/2. 

The Minkowski sausage In is a polygon in 1? defined as 
follows: (i) Iq is a square of side one. (ii) I„+i is obtained 
by replacing each side of I„ by a path of length three. The 



angles in the path are determined by the position of the side 
in In- The first and third segments of the path follow the 
direction of the replaced side. The two angles are first left 
then right. Analogously, we can construct In+i by attaching 
to the sides of I„ four of its copies (see Fig.|2l). The polygon 
In can be used to tessellate the plane. From the definition, we 
can determine p (I„) =4x3" and a (!„) = 2 x 3" — 2. Here 
too we have S{n) — p{n) — a{n). Hence, 7 (!„) = 1/2. 

The Moore polygon is a "closed version" of the Moore 
curve. It is a polygon in I? defined by a closed path expressed 
as an L-system. A Lindenmayer system (for short, L-system) 
H is a quadruple {V, C, A, R), where V is a set of variables, 
C a set of constants, A a set of axioms, and B a set of pro- 
duction rules. An L-system allows the recursive construction 
of words (or, equivalently, sequence of symbols) whose let- 
ters are elements from V and C. An axiom is a word at time 
t = 0. At each time step t+1, the production rules are applied 
to the word given by the L-system at time t. Only variables 
are replaced according to the production rules. On the basis 
of these definitions, we can write Ain = {V, C, A, R), where 
V = {a, b}, C = {+, -}, A = {aFa + F + aFa}, and 
R= {a^ -bF + aFa + Fb-, b +aF - bFb - Fa+}. 
The letter F indicates a segment of length one in I?. The first 
segment of A4o specified by the axiom in A is {(0, 0), (1, 0)}. 
The symbols + and — mean "turn left in Z^" and "turn right in 
Z^", respectively. The sequences — h and H — have no mean- 
ing and can be deleted. For instance, the polygon A^i is then 
given by the the following word: —bF+aFa+Fb—F—bF+ 
aFa + FbFbF + aFa +Fb-F-bF + aFa + Fb~F. Notice 
that in order to close 1 we need to replace ■ ■ ■ + Fb — F with 
• • • + FbF in the obtained word. This operation is required 
for every n. Once we have generated the polygon, we blow 
it up by replacing each square of side one with a square com- 
prising four of its copies. The occurrences of letter F in the 
word produced by 1 is 16. In general, the number of occur- 
rences of F in the word produced by Ain equals the perime- 
ter of Mn- From the definition, this is p{Mn) = 2 x 4"+^, 
taking into account the blowing up operation. The number 
of — ("turn right") symbols, excluding the initial one, in the 
word produced by A^„, is exactly equal to the number of in- 
ward angles of A4„: a{Mn) = | (-1)" + §4" - 2. From 
S — p{Mn) — a{Mn), we can compute 7 (Mn) = 4/5. 

The Vicsek snowflake on 1? , denoted by V„, is a bounded 
region of 1? defined iteratively as follows: {i) Vg is a single 
1x1 square, {ii) We obtain V„+i by attaching 4 copies of 
V„ to its corners (see Fig. |2]i. Each square comprising V„ 
has side one. For this fractal we have p (Vn) = 20 x 5"^^ 
and a (V„) = 2 x 5" — 2. The number of adjacent squares is 
S{n) = p{n) — a{n), which gives 7(V) = ^■ 

The quadratic Koch polygon, /C„, is a polygon in Z^ based 
on the Koch curve. Essentially, it consists of a region bounded 
by two mirroring copies of the Koch curve. As the Moore 
polygon, ICn is defined by an L-system and specified by a 
path. The path giving rise to ICq is given axiomatically as 
{(0, 0), (1, 0)}. Then ICq is a square of side one. The produc- 
tion rule is F^F + F — F — F + F, where F indicates 
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TABLE I: Fractal entanglement 7, perimeter p{n), entropy of entan- 
glement S{n) for a state in £■ for several fractal bipartitions {A, B) 
of the square lattice. Here D is the Hausdorff dimension of the curve 
separating the regions An and B„. Forp(n) and S{n) only the lead- 
ing term is shown. 



Fractal 


1 






D 


1. Sierpinski carpet 


99 
224 


5^ 


99 on 
280° 


log 8 
log 3 


2. Greek Cross 


1 
2 


8n 


4n 


2 


3. Minkowski Sausage 


1 
2 


4x3" 


2x3" 


logs 
log 3 


4. Vicsek Snowflake 


1 
2 


4x5" 


2x5" 


log S 
log 3 


5. Quadratic Koch 


58 
125 


4x5" 


232 rn 
125° 


log S 
log 3 


6. Moore Polygon 


4 
5 


2 x 4"+i 


32 An 
5 


log 9 
log 6 


7. T-Square 


1 
2 


16 X 3" 


92 on 
9 ^ 


2 


8. Chessboard 


1 
4 


8n^ 


2n^ 


2 



again a segment of length one in 1? . The fractal has a pat- 
tern similar to that of the Vicsek snowflake and indeed has the 
same Hausdorff dimension (see Table Nevertheless, the 
results for the scaling of the entanglement are different. The 
perimeter can be computed as p{n) = 4 x 5". The number 
h{n) of holes is /i = x 5" + ^3" - 1, for n > 3. One 
can easily see that a = (p — Ah)/ 2 and therefore from Eq. ^ 
S = ^-h = f§5" - i3" + 1. In the Hmit of large n, we 
obtain 7 = 58/125. 

The T-square polygon on 1?, £„, is obtained by superim- 
posing four copies of f„_i on the corners of a square of side 
2n+i -pjjg covered by each copy is exactly a square of 
side 2". The perimeter of is p (£„) = 16 x 3" - 8 x 2". We 
have S{£q) = 4, S {£1) = 24, and 5 (£„) = 3S'(f„_i) + 
2"+i-8 = f 3"+2"+i-8+24xS'(n, 3) = f 3"-4x2"+4, 
where 3) := (1 + 3""^ _ 2"-i)/2 is the n-th Stirling 
number of the second kind. Hence, 7 = 1/2. 

The chessboard C„ is the bounded region of 1? defined as 
follows. Let Ci be a 2 X 2 square with two holes in the upper 
right and bottom left corner Then C„+i is obtained by placing 
4 copies of Cn on all the quadrants of a 2" x 2" square on 1?. 
The perimeter is p = 2n. The number of adjacent squares is 
exactly h = n/2. Therefore it is immediate that 7 = Ns/p = 
1 /4 for every size n. It is obvious that this is a lower bound for 
the entanglement on the square lattice for a state in C, since 
the chessboard maximizes the number of holes of side 1. 

Conclusions. — This work has, for the first time, explored 
the relationship between entanglement entropy and the frac- 
tality of the bipartition in a spin system. We have calculated 
the scaling of entanglement S with the length p of the bound- 
ary in the ground state of the Z2 topological phase associated 
with the toric code, for various fractal boundaries. We have 
shown that this provides an upper bound on the entanglement 
in the entire topological phase. Unlike the case of a regular 
boundary, the ratio 7 — S/p for large p is not exactly 1 but a 
smaller fraction, so that the general bound for the area law is 
still obeyed. The fractal nature of the bipartition is revealed in 
the total amount of entanglement present in the system. There 
is less entanglement in a fractal bipartition. We also found 



that the ratio 7 is always at most the inverse of the Hausdorff 
dimension D. We conjecture this last claim to hold in general 
for topologically ordered states. Moreover, different fractals 
with the same Hausdorff dimension can have different 7, so 
that this is a useful quantity to classify fractals with. We chose 
the toric code because in this case it is simple to compute the 
entanglement. It would be interesting to consider other types 
of topologically ordered states and explore whether the behav- 
ior we have observed is general for any quantum system with 
finite correlation length. Finally, since the scaling of entan- 
glement with the boundary of the system is less than 1, we be- 
lieve that a renormalization group algorithm based on blocks 
of spins that grow like fractals, might be potentially more ef- 
ficient. Indeed, in this regard the chessboard appears to be the 
most attractive of all the fractals we have considered. 
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